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In the leading twist approximation of the Wilson operator product expansion 
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with "frozen" and analytic strong coupling constants we show that Bessel-inspired 
behavior of the structure function F2 at small x, obtained for a flat initial condition 



' i n the DGLAP evolution equations, leads to good agreement with the deep inelastic 

ft, , scattering experimental data from HERA. 
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. The experimental data from HERA on the deep-inelastic scattering (DIS) structure 

tion (SF) F 2 fl, and its derivatives dF 2 /d\n(Q 2 ) Q, Q and d\nF 2 /d\n(l/x) fl, 2, 3 
O ■ 

enter us into a very interesting kinematical range for testing the theoretical ideas on the 
behavior of quarks and gluons carrying a very low fraction of momentum of the proton, the 
so-called small- 2 region. In this limit one expects that the conventional treatment based on 
the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equations does not account for 
contributions to the cross section which are leading in a s ln(l/x) and, moreover, the parton 
densities (PD), in particular the gluons, are becoming large and need to develop a high 
density formulation of QCD. 

However, the reasonable agreement between HERA data and the next-to-leading-order 
(NLO) approximation of perturbative QCD has been observed for Q 2 > 2 GeV 2 (see reviews 
in 6| and references therein) and, thus, perturbative QCD could describe the evolution of 
F 2 and its derivatives up to very low Q 2 values, traditionally explained by soft processes. 
The standard program to study the x behavior of quarks and gluons is carried out by 



comparison of data with the numerical solution of the DGLAP equations [7.] 1 by fitting 
the parameters of the x-profile of partons at some initial Q$ and the QCD energy scale A 



tnepE 

0,Q, 



However, for analyzing exclusively the small-x region, there is the alternative of 



doing a simpler analysis by using some of the existing analytical solutions of DGLAP in the 



small-x limit [3]-[l5j. This was done so in 12j where it was pointed out that the HERA 
small-x data can be interpreted in terms of the so-called doubled asymptotic scaling (DAS) 
phenomenon related to the asymptotic behavior of the DGLAP evolution discovered many 
years ago flf] ]. 

The study of [l2] was extended in Q, Q to include the finite parts of anomalous 



dimensions of Wilson operators and Wilson coefficients 2 . This has led to predictions [jjj] 
of the small-x asymptotic PD form in the framework of the DGLAP dynamics starting at 
some Qq with the flat function 



f a {Ql) = K (hereafter a = q, g), 



(1) 



where f a are the parton distributions multiplied by x and A a are unknown parameters to 
be determined from the data. 
We refer to the approach of 
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151 ] as generalized DAS approximation. In the 



approach the flat initial conditions in Eq. (TjQ) determine the basic role of the singular 
parts of anomalous dimensions, as in the standard DAS case, while the contribution from 
finite parts of anomalous dimensions and from Wilson coefficients can be considered as 
corrections which are, however, imp ortant for better agreement with experimental data. In 



the present paper, similary to 
quark component. 



121- 



15], we neglect the contribution from the non-singlet 
is supported by the actual experi- 



The use of the flat initial condition given in Eq. 

nnnn 

mental situation: low-Q 2 data |l|, |3|, ll7|, [18] are well described for Q 2 < 0.4 GeV 2 by Regge 



theory with Pomeron intercept «p(C 
one. The small rise of HERA data 



\p + 1 = 1.08, closed to the standard («p(0) = 1) 



18 



19] at low Q 2 can be explained, for example, 



by contributions of higher twist operators (see 15|). 

The purpose of this paper is to demostrate a good agreement between the predictions 



1 At small x there is another approach based on the Balitsky-Fadin-Kuraev-Lipatov (BFKL) equation 

of this work. 



whose application is out of the scope 
In the standard DAS approximation 



16| only the singular parts of the anomalous dimensions were used. 



from the generalized DAS approach and the HERA experimental data 



'see Figs. 1-3) 



for SF F 2 and a 
and ZEUS data 



so to compare the predictions for the slope 91nF 2 /91n(l/x) with the HI 



sl, Q, S]. Looking at the HI data [l] points shown in Figs. 4 and 5 one can 
conclude that \{Q 2 ) is independent on x within the experimental uncertainties for fixed Q 2 
in the range x < 0.01. Indeed, the data are well described by the power behavior 



F 2 (x,Q 2 ) = Cx~^\ 



(2) 



where X(Q 2 ) = aln(Q 2 /A 2 ) with C « 0.18, a « 0.048 and A = 292 MeV [4|. The linear rise 
of the exponent \{Q 2 ) with InQ 2 is also explicitly shown in Figs. 4 and 5 by the dashed 
line. 

The rise of \{Q 2 ) linearly with \\iQ 2 can be tracted in strong nonperturbative way (see 
[20I ] and references therein), i.e., \{Q 2 ) ~ l/a s (Q 2 ). The previous analysis 21], however, 
demonstrated that the rise can be explained naturally in the framework of perturbative 
QCD. 

The ZEUS and HI Collaborations have also presented new preliminary data for 

X(Q 2 ) at quite low values of Q 2 . As it is possible to see in Fig. 8 of jsj], the ZEUS value 
for X(Q 2 ) is consistent with a constant ~ 0.1 at Q 2 < 0.6 GeV 2 , as it is expected under the 
assumption of single soft Pomeron exchange within the framework of Regge phenomenology. 
These points lie slightly below the corresponding ZEUS data but all the results are in 
agreement within modern experimental errors. 

It is important to extend the analysis of [21I to low Q 2 range with a help of well-known 
infrared modifications of the strong coupling constant. We will use the "frozen" and analytic 
versions (see, [22J and [23(, respectively). 



The paper is organized as follows. Section 2 contains basic formulae, which are needed 
for the present study and were previously obtained in 21|. In Section 3 we compare 



our calculations with HI and ZEUS experimental data and present the obtained results. 
Some discussions can be found in Conclusions. 



II. GENERALIZED DAS APPROACH 



The flat initial condition ([T]) corresponds to the case when parton distributions tend to 



4 



some constant value at x — > and at some initial value Qq. The main ingredients of the 
results HQ , are: 

• Both, the gluon and quark singlet densities are presented in terms of two components 
(" + " and " — " ) which are obtained from the analytic Q 2 -dependent expressions of 
the corresponding (" + " and " — " ) PD moments. 

• The twist-two part of the " — " component is constant at small x at any values of Q 2 , 
whereas the one of the " + " component grows at Q 2 > Qq as 

- e CT , a = 2 



d+s - lD+ + d + j-\ p 



where a and p are the generalized Ball-Forte variables, 

s=ln feiD' '=«-«*>-«-w>- D ^w (4) 

Hereafter we use the notation a s = a s j (47r). The first two coefficients of the QCD /3-function 
in the MS-scheme are (3q = 11 — (2/3)/ and f3\ = 102 — (114/9)/ with / is being the number 
of active quark flavors. 

Note here that the perturbative coupling constant a s (Q 2 ) is different at the leading-order 
(LO) and NLO approximations. Indeed, from the renormalization group equation we can 
obtain the following equations for the coupling constant 

1 a ,J Q 2 X 

at the LO approximation and 



1 (3 : 



a s (Q 2 ) A 



i3 + f3ia s (Q 2 ) 



2' 



( 6 ) 



at the NLO approximation. Usually at the NLO level MS-scheme is used, so we apply 



A = Ajyfg below. 



A. Parton distributions and the structure function F2 

Here, for simplicity we consider only the LO approximation 3 . The results for parton 
densities and F 2 are following: 



The NLO results may be found in 
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The structure function F 2 has the form 

l2\ „ n f rfZ 



F 2 (x,Q 2 ) = ef q (x,Q%, 
f a (x,Q 2 ) = f^(x,Q 2 ) + f-(x,Q 2 ), (7) 

where e = (J2i e l)/ f is the average charge square. 



• The small-x asymptotic results for PD are 



f+(x,Q 2 ) = (A g + ^A q y o (a)e-~ d +U s + O(p), (8) 



ft^Q 2 ) = {{A g + A -A q )p~h{a)e- d ^ + O(p), (9) 
4 



f g (x,Q 2 ) = --V" d - (1)s + 0(x), (10) 
f~(x,Q 2 ) = A g e- d '^ s + 0(x), (11) 

where d + (l) = 1 + 20//(27/? ) and d_(l) = 16//(27/3 ) are the regular parts of the 
anomalous dimensions d + (n) and d_(n), respectively, in the limit n — > l 4 . Here n is 
the variable in Mellin space. The functions l v [y = 0, 1) are related to the modified 
Bessel function I v and to the Bessel function J v by: 

{ IJa), if s > 

W = " • (12) 

[ i-"J u (ia), i 2 = -1, if s < 

At the LO, the variables a and p are given by Eq. ([3]) when p = 0. 



B. Effective slopes 

Contrary to the approach in Jj], 
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15| various groups have been able to fit the 
available data using a hard input at small x: x _A , A > with different A values at low and 
high Q 2 (see 



24 



25 



26 



27 



28 



29j). Such results are well-known at low Q 2 values 25] 



At large Q 2 values, for the modern HERA data it is also not very surprising, because they 
cannot distinguish between the behavior based on a steep input parton parameterization, at 
quite large Q 2 , and the steep form acquired after the dynamical evolution from a flat initial 
condition at quite low Q 2 values. 



4 We denote the singular and regular parts of a given quantity fc(n) in the limit n — > 1 by k(n) and k(n), 
respectively. 



As it has been mentioned above and shown in Il5| . the behavior of parton densities 
and F 2 given in the Bessel-like form by generalized DAS approach can mimic a power law 
shape over a limited region of x and Q 2 



f a (x,Q 2 ) ~ x~ x ° S( - x ' Q2) and F 2 (x,Q 2 ) ~ x~ X ^ 



The effective slopes \ e a (x, Q ) and A|P (x, Q ) have the form: 
X°*(x Q 2 ) = f ^ Q2) ,M 

Aefffx = / g + (^ 2 ) / 2 (a)(l ~ 20a,(g 2 )) + 20a s (Q 2 )I 1 (a)/p 
9 ) f q (x,Q 2 ) P h(a)(l-2Va s (Q 2 )) + 2Va s (Q 2 )Ua)/p 

, eS( n2) Ag"(s, Q 2 ) fffo Q 2 ) + (2/)/3as(g 2 ) Ag"(s, Q 2 ) Q 2 ) 
F2lX ' V J /,(*, Q 2 ) + (2/)/3a s (Q 2 ) g 2 ) ' 



where the exact form of parton densities can be found in 

The results (1T31) (and also (|T4l) -(fT6l ) below) are given at the NLO approximation. To 
obtain the LO one, it is necessary to cancel the term ~ a s (Q 2 ) and to use Eqs. (l8l)-( |TTi) for 
parton densities f a (x,Q 2 ). 

The effective slopes A„ ff and Af? depend on the magnitudes A a of the initial PD and 
also on the chosen input values of Ql and A. To compare with the experimental data it 
is necessary the exact expressions ( fl3l) . but for qualitative analysis it is better to use an 
approximation. 



C. Asymptotic form of the effective slopes 

At quite large values of Q 2 , where the "— " component is negligible, the dependence on 
the initial PD disappears, having in this case for the asymptotic behavior the following 
expressions 5 : 

Af ^ g2) ^|M Kp __i_, (14) 

5 The asymptotic formulae given in Eqs. (TH]) ~ (ri!))) work quite well at any Q 2 > Qq values, because at 
Q 2 = Qq the values of A® ff and are equal zero. The use of approximations in Eqs. (jT3J)— dXHl) instead 
of the exact results given in Eq. (fiU)) underestimates (overestimates) only slightly the gluon (quark) slope 
at Q 2 > Ql 
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v e ff ,as fx Q2) M°)(l - 20a 5 (g 2 )) + 20a a (Q*)I 1 (<T)/p 
q ^ ) P h(a)(l-20a s (Q2)) + 20a s (Q*)Io(a)/p 



3 | 1Qa °(Q 2 ) (15) 

P 41n(l/x) pln(l/x)' 1 ' 

A£-(*, Q 2 ) = P ! M + 2Qa s (Q 2 ) ( 1 - toWl'W ) 

41n(l/x) pm(l/x) y pm(l/x) 
where the symbol ~ marks the approximation obtained in the expansion of the usual and 
modified Bessel functions in (I12j) . These approximations are accurate only at very large a 
values (i.e. at very large Q 2 and/or very small x). 

As one can see from Eqs. (fT4j) and (jT5]h the gluon effective slope A^ is larger than the 



quark slope A^ ff , which is in excellent agreement with MRS 30j and GRV 10j analyses (see 

also 

We would like to note that at the NLO approximation the slope Ap^' a& (x, Q 2 ) lies between 
quark and gluon ones but closely to quark slope A^ ff ' as (x, Q 2 ). Indeed, 



A«'»(x,Q') - Af.»(x,Q 3 ) « -j^. (18) 

Both slopes A^ ff (x,Q 2 ) decrease with increasing x (see Fig. 5). A x-dependence of 
the slope should not appear for a PD within a Regge type asymptotic (x~ A ) and precise 
measurement of the slope X^(x, Q 2 ) may lead to the possibility to verify the type of small- x 
asymptotics of parton distributions. 

III. COMPARISON WITH EXPERIMENTAL DATA 



Using the results of previous section we have analyzed HERA data for Fo and the slope 
<91nF 2 /<91n(l/x) at small x from the HI and ZEUS Collaborations [l, 2, ^,0,0]. 

In order to keep the analysis as simple as possible, we fix / = 4 and a s (M§) = 0.1166 
(i.e., A^ 4 ) = 284 MeV) in agreement with the more recent ZEUS results {2J. 

As it is possible to see in Figs. 1-3 (see also Q, Q), the twist-two approximation is 
reasonable at Q 2 > 2 GeV 2 . At smaller Q 2 , some modification of the approximation should 
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be considered. In the recent article 15j we have added the higher twist corrections. For 
renormalon model of higher twists, we have found a good agreement with experimental data 
at essentially lower Q 2 values: Q 2 > 0.5 GeV 2 (see Figs. 2 and 3). 

Moreover, the results of fits in 15] have an important property: they are very similar in 
LO and NLO approximations of perturbation theory. The similarity is related to the fact 
that the small-x asymptotics of the NLO corrections are usually large and negative (see, for 



example, a s -corrections 3l| to BFKL kernel [8j 6 ). Then, the LO form ~ ce s (Q 2 ) for some 
observable and the NLO one ~ a s (Q 2 )(l — Ka s (Q 2 )) with a large value of K are similar, 
because A 3> Alo 7 and, thus, a s {Q 2 ) at LO is considerably smaller then a s {Q 2 ) at NLO for 
HERA Q 2 values. 

In other words, performing some resummation procedure (such as Grunberg's effective- 
charge method [32]), one can see that the NLO form may be represented as ~ a s (Q 2 s 
where Q 2 S >■ Q 2 . Indeed, from different studies 



33 



34 



35] , it is well known that at small-x 
values the effective argument of the coupling constant is higher then Q 2 . 

Here, to improve the agreement at small Q 2 values, we modify the QCD coupling constant. 
We consider two modifications, which effectively increase the argument of the coupling 



34 



35|) 



constant at small Q 2 values (in agreement with [3 

In one case, which is more phenomenological, we introduce freezing of the coupling con 
stant by changing its argument Q 2 — > Q 2 + M 2 , where M p is the p- meson mass (see 22j ) 
Thus, in the formulae of the Section 2 we should do the following replacement: 



a s (Q 2 ) - a fr (Q 2 ) = a s (Q 2 + M 2 



(19) 



The second possibility incorporates the Shirkov-Solovtsov idea [23j, [36|, |37[ about ana- 
lyticity of the coupling constant that leads to the additional its power dependence. Then, 
in the formulae of the previous section the coupling constant a s (Q 2 ) should be replaced as 
follows: 



a^(Q 2 ) = a s (Q 2 ) - 



1 



A 2 

A LO 



A) Q 2 - Aio 



(20) 



6 It seems that it is a property of any processes in which gluons, but not quarks play a basic role. 

7 The equality of a s (M^) at LO and NLO approximations, where Mz is the Z-boson mass, relates A and 
A LO : AM = 284 MeV (as in [2]) corresponds to A LO = 112 MeV (see 0). 



9 



at the LO approximation and 

= (21) 

at the NLO approximation, where the symbol . . . stands for terms which have negligible 



23) 



contributions at Q > 1 GeV 

Figure 4 shows the experimental data for X e p 2 (x,Q 2 ) 9 at x ~ 10~ 3 , which represents an 
average of the x-values of HERA experimental data. The top dashed line represents the 
aforementioned linear rise of X(Q 2 ) with ln(Q 2 ). 

Figure 4 demonstrates that the theoretical description of the small-Q 2 ZEUS data for 
X c p 2 (x, Q 2 ) by NLO QCD is significantly improved by implementing the "frozen" and analytic 
coupling constants af T (Q 2 ) and a an (<5 2 ), respectively, which in turn lead to very close results 



(see also |39l|). 

Indeed, the fits for F 2 {x,Q 2 ) in [l5j yielded Q 2 « 0.5-0.8 GeV 2 . So, initially we had 
Xf 2 (x,Ql) = 0, as suggested by Eq. (□). The replacements of Eqs. (TT9|) . (1201) and (I2T1) 
modify the value of X e p 2 (x, Qq). For the "frozen" and analytic coupling constants a^{Q 2 ) and 
a an (Q 2 ), the value of X c p 2 (x, Qq) is nonzero and the slopes are quite close to the experimental 
data at Q 2 « 0.5 GeV 2 . Nevertheless, for Q 2 < 0.5 GeV 2 , there is still some disagreement 
with the data, which needs additional investigation. 

Figure 5 shows the x-dependence of the slope X e p 2 (x,Q 2 ). One observes good agreement 
between the experimental data and the generalized DAS approach for a broad range of 
small-x values. The absence of a variation with x of Af?(x, Q 2 ) at small Q 2 values is related 
to the small values of the variable p there. 

At large Q 2 values, the x-dependence of X e p 2 (x, Q 2 ) is rather strong. However, it is well 
known that the boundaries and mean values of the experimental x ranges [4| increase propor- 
tionally with Q 2 , which is related to the kinematical restrictions in the HERA experiments: 
x ~ 10~ 4 x Q 2 (see , and, for example, Fig. 1 of [sj]). We show only the case with 
the "frozen" coupling constant because at large Q 2 values all results are very similar. 



Note that in 



36l. l37| more accurate, but essentially more cumbersome approximations of a an (Q ) have 
been proposed. We limit ourselves by above simple form (|20|) . (f2Tj) and plan to add the other modifications 
in our future investigations. 
9 Using the "frozen" and analytic coupling constants, the experimental data for F<z(x, Q 2 ) have been anal- 



ysed recenly in [38 1. 
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From Fig. 5, one can see that HERA experimental data are close to X e p 2 (x,Q 2 ) at x ~ 
10~ 4 -=- 10~ 5 for Q 2 = 4 GeV 2 and at x ~ 10~ 2 for Q 2 = 100 GeV 2 . Indeed, the correlations 
between x and Q 2 in the form x c s = a x 10~ 4 x Q 2 with a = 0.1 and 1 lead to a modification 



of the Q 2 evolution which starts to resemble \nQ 2 , rather than lnlnQ 2 as is standard 21]. 



IV. CONCLUSIONS 



We shown the Q 2 -dependence of the structure function F 2 and the slope Af? = 
d hi F 2 / d ln(l / x) at small-x values in the framework of perturbative QCD. Our twist-two 
results are in very good agreement with precise HERA data at Q 2 > 2 GeV 2 , where per- 
turbative theory can be applicable. The application of the "frozen" and analytic coupling 



SB, 



constants af T (Q ) and a &n {Q ) improves the agreement with the recent HERA data 
for the slope X c p 2 (x,Q 2 ) for small Q 2 values, Q 2 > 0.5 GeV 2 . 

As a next step of investigations, we plan to fit the HERA data l|, Q, of F 2 {x,Q 2 ) 

directly, using the "frozen" and analytic coupling constants in both the LO and NLO 
approximations, in order to improve the agreement with HERA data at small Q 2 values. 
Several versions of the analytical coupling constant will be used. 

This work was supported by RFBR grant 07-02-01046-a. A.V.K. thanks the Organizing 
Committee of Workshop on Physics of Fundamental Interactions for invitation. 
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Figure captions 



Fig. 1. Fo(x,Q ) as a function of x 



are from HI [1] (open points) and ZEUS 2] (solid points) at Q 2 > 1.5 GeV 2 . The solid 
curve represents the NLO fit. The dashed curve (hardly distinguishable from the solid one) 
represents the LO fit. 

Fig. 2. F2(x,Q 2 ) as a function of x for different Q 2 bins. The experimental points are 
same as on Fig. 1. The solid curve represents the NLO fit. The dash-dotted curve represents 
the BFKL-motivated estimation for higher-twist corrections to F2(x,Q 2 ) (see [151]). The 
dashed curve is obtained from the fits at the NLO, when the renormalon contributions of 
higher-twist terms have been incorporated. 

Fig. 3. F 2 (x,(5 2 ) as a function of x for different Q 2 bins. The experimental points 
are from HI [1] (open points) and ZEUS 2] (solid points) at Q 2 > 0.5 GeV 2 . The solid 
curve represents the NLO fit. The dashed curve is from the fits at the NLO with the 
renormalon contributions of higher-twist terms incorporated. The dash-dotted curve (hardly 
distinguishable from the dashed one) represents the LO fit with the renormalon contributions 
of higher-twist terms incorporated. 

Fig. 4. The values of effective slope Af? as a function of Q 2 for x = 10 -3 . The 
experimental points are from HI 4, 5] (open points) and ZEUS jjj (solid points). The 
solid curve represents the NLO fit. The dash-dotted and lower dashed curves represent the 
NLO fits with "frozen" and analytic coupling constants, respectively. The top dashed line 
represents the fit from 4|. 

Figure 5. The values of effective slope Af? as a function of Q 2 . The experimental points 
are same as on Fig. 4. The dashed line represents the fit from [4]. The solid curves represent 
the NLO fits with "frozen" coupling constant at x — 10~ 2 and x = 10~ 5 . 
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FIG. 1: F2(x, Q 2 ) as a function of x for different Q 2 bins. The experimental points are from HI [l| 
(open points) and ZEUS [2] (solid points) at Q 2 > 1.5 GeV 2 . The solid curve represents the NLO 
fit. The dashed curve (hardly distinguishable from the solid one) represents the LO fit. 
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FIG. 2: F2(x,Q 2 ) as a function of x for different Q 2 bins. The experimental points are same 
as on Fig. 1. The solid curve represents the NLO fit. The dash-dotted curve represents the 
BFKL-motivated estimation for higher-twist corrections to F2(x,Q 2 ) (see [15l]). The dashed curve 
is obtained from the fits at the NLO, when the renormalon contributions of higher- twist terms 
have been incorporated. 




G. 3: F2(x,Q 2 ) as a function of x for different Q 2 bins. The experimental points are from HI 
l[ (open points) and ZEUS [2] (solid points) at Q 2 > 0.5 GeV 2 . The solid curve represents the 
NLO fit. The dashed curve is from the fits at the NLO with the renormalon contributions of 
higher-twist terms incorporated. The dash-dotted curve (hardly distinguishable from the dashed 
one) represents the LO fit with the renormalon contributions of higher-twist terms incorporated. 
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Q 2 



FIG. 4: The values of effective slope A|? as a function of Q 2 for x = 10 -3 . The experimental points 
are from HI 4, 5] (open points) and ZEUS j| (solid points). The solid curve represents the NLO 
fit. The dash-dotted and lower dashed curves represent the NLO fits with "frozen" and analytic 
coupling constants, respectively. The top dashed line represents the fit from (4]. 
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FIG. 5: The values of effective slope as a function of Q 2 . The experimental points are same 
as on Fig. 4. The dashed line represents the fit from [4|. The solid curves represent the NLO fits 
with "frozen" coupling constant at x = 10 -2 and x = 10~ 5 . 



